Jamming of frictional tetrahedra 
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We prepare packings of frictional tetrahedra with volume fractions <j) ranging from 0.470 to 0.612 
using three different experimental protocols under isobaric conditions. Analysis via X-ray micro- 
tomography reveals that the contact number Z grows with <j), but does depend on the preparation 
protocol. While Z saturates at high packing fraction, the number of constraint per particle C in- 
creases further, supported by an increasing number of edge-to-face and face-to-face contacts. Even 
the loosest packings are strongly hyperstatic i.e. mechanically over-determined with C approximately 
twice the degrees of freedom each particle possesses. 

PACS numbers: 83.80.Fg, 45.70.Cc 




FIG. 1. From left to right: Close-up view of the surface 
of an experimental tetrahedra packing; Rendering of X-Ray 
Computed- Tomography data; Detected tetrahedra, colored 
by orientation. (Color online) 



Introduction. - The jamming paradigm [Tl[2] provides 
an unified viewpoint on the onset of mechanical stability 
in particulate soft matter systems like wet foams, emul- 
sion, and colloids. For such frictionless spherical par- 
ticles there exists a critical volume fraction 0j 0.64 
where the system becomes isostatic, i.e. the number of 
constraints imposed by the contacts corresponds to the 
degrees of freedom of the particles. Higher densities can 
be achieved by increasing the pressure which leads to ex- 
cess contacts; these excess contacts then determine the 
mechanical properties of the system. 

Granular matter differs in two ways from these sys- 
tems. First, its nonzero surface friction allows for tan- 
gential forces at particle contacts. These result in me- 
chanical stable packings over a full range of </> at any 
given pressure. Such packings differ in their structure, 
not their compression. However, one can still define a 
meaningful distance to isostaticity [3]. 

Secondly, most real world granular materials like sand, 
sugar or soil are not spherical. Only quite recently the 
jamming approach been been applied to frictionless ellip- 
soids PI^S] and tetrahedra [B]. Tetrahedra provide besides 



jamming a number of other interesting packing problems 
like the existence of quasi crystals [7] or the quest for 
the densest packing [5HTU] , which presently is believed to 
be a dimer crystal with </) « 0.8563 nj[]. However, ex- 
perimental packings of necessarily frictional tetrahedra 
explore a much lower range of (f) [121 . 

From a jamming point of view tetrahedra introduce 
a new element: unlike spheres they have four different 
types of contacts: face to face, edge to face, edge to edge, 
and vertex to face. Nevertheless, previous numerical re- 
sults have found frictionless uncompressed packings of 
tetrahedra to be isostatic [B1[T21[I5; the same result has 
been reported by the only experiment measuring contact 
numbers on frictional tetrahedral dice |15j . 

In this letter we use X-ray tomography to determine 
the particle positions and orientations inside tetrahedra 
packings. We show that even our loosest packings are 
highly hyperstatic and that the contact number are in 
general protocol-dependent. 

Experimental setup. - The tetrahedra particles in our 
study are produced by injection moulding of polypropy- 
lene, and have a side length of a = 7 mm (see Fig. la I. 
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The typical radius of curvature of edges and corners is 
150 ± 50 /im which corresponds to 2.5% of a. All parti- 
cles have small mold marks on one of the 6 edges, the 
height of this protrusions amounts to 300 /im at max- 
imum [TCj. The particles have a coefficient of friction 
jjL = 0.87 ± 0.03, which is measured using a tilted plane 
and finding the angle at which face-to-face contacts start 
to slide past each other. 

In order to explore the structure of tetrahedral pack- 
ings at different densities, we compactify them by verti- 
cal shaking. We start from initially loose packings which 
are prepared by filling 14406 particles via a slowly lifted 
funnel into a cylinder of inner diameter D = 104 mm 
(w 15 a). This method creates packings with (j) between 
0.470 and 0.477, the lowest values found in this study. 
An electromagnetic shaker (LDS V555) then applies sine 
pulses ("taps") with a width of 50ms and a repetition 
rate of 3 Hz. The peak acceleration Ap per tap is mea- 
sured at the base plate of the container (Kistler 8763B) 
and reported in units of the dimensionless acceleration 
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r = Ap/g (where g = Q.Slm/s^). 

A laser distance scanner (MicroEpsilon ILD1402) 
mounted on a translation stage measures the surface 
height profile after shaking; from this profile the corre- 
sponding (j) is computed using calibration measurements 
from tomographic reconstruction. This results in an 
accuracy of the reported values of ± 0.003. Tomo- 
grams of the prepared packings are acquired using X-ray 
Computed Tomography (GE Nanotom) with a tungsten 
target and 120 kV acceleration voltage. Each dataset 
comprises 1152x1152x1076 voxels with a resolution of 
100 /j,m/voxel. 

Preparation protocols. - We use three different prepa- 
ration protocols. The first one will be referred to as TAP: 
Starting from a loose packing, a number of taps ranging 
from 20 to 6400 with a constant acceleration F = 4 is 
applied. As Fig. [2a| shows, the packing fraction increases 
monotonely with the number of taps. Fig. |2b| demon- 
strates that the system reaches after a sufficient number 
of taps a history-independent steady state with a packing 
fraction <j)ss- 

Our second protocol RAMP is intended to study the 
possible range of steady state volume fractions ^s^. Fol- 
lowing Nowak et a/.|17j. the tap intensity F is increased 
stepwise from 0.5 to 10, decreased back to 0.5, and in- 
creased to 10 again. At each step, 10* taps are applied 
and the resulting <f>ss is measured. The evolution of (p^s is 
depicted in Fig. [2cj and comprises the following regimes: 
The loose packing irreversibly compacts when F is in- 
creased from 0.5 to 3, but dilates again when the ac- 
celeration is further ramped to up F = 10, leading to 
4>ss ~ 0.55. Subsequently, decreasing the tap intensity 
down to F = 0.5 creates a dense packing with 0ss ~ 0.60. 
Increasing F again to 10.0 demonstrates the reversibil- 
ity i.e. no hysteresis. We refer to the end points of the 
reversible branch as RAMPr=o.5 and RAMPr=io in our 
tomograhic analysis. 

The third preparation protocol VIB is inspired by epi- 
taxial growth: particles are deposited at a rate of approx. 
15/second, while continuous sine vibration at 100 Hz and 
F = 10 is applied; a similar method has been employed 
to create dense packings of ellipsoids [H]. Protocol VIB 
creates the densest packings in our study with — 0.612. 
Apparently the slow deposition and the continuous agi- 
tation of the surface layer enables the particles to explore 
more effectively the configurational phase space and thus 
find configurations with lower potential energy. 

Our range of mechanically stable packings of </) = 0.470 
to 0.612 is in accordance with the range of = 0.48 to 
0.64 found by Baker and KudroUi ^ for tetrahedral dice 
and ceramic tetrahedra particles prepared by shaking and 
a water-fluidized bed. Our volume fractions are however 
significantly lower than the value of (j) ~ 0.76 reported 
by Jaoshvili et al. [TS] , who prepared packings of tetrahe- 
dral dice under water with an unspecified protocol. This 
discrepancy is unlikely due to the larger corner radius of 
the dice as optimal packings of weakly truncated tetrahe- 
dra (resembling the dice) do not pack denser than ideal 
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FIG. 2. (a) Protocol TAP: Tapping at F = 4 compacts the 
packing, (b) Packing fraction evolution of an initially loose 
vs. an initially dense packing, (c) Protocol RAMP: At each 
acceleration step, the packing is tapped 10* times. Open cir- 
cles (green): initial, irreversible compaction. Open downward 
triangles (orange): decreasing acceleration. Filled upward tri- 
angles (violet): increasing acceleration. Packings at points 
RAMPr=o.5 and RAMPr=io are analyzed via tomography. 
All errorbars are averaged over 4 independent runs. (Color 
online) 



tetrahedra|19|. Judging from our experience with spheres 
also neither the different values of /i, nor the presence of 
interstitial water (c.f. |20] for a demonstration of static 
friction under water) can explain this large difference. 

Tomography analysis. - Tomographic reconstruction 
of the packings allows us to investigate the geometrical 
properties of the packing by registering the positions and 
orientations of all tetrahedra. The numerical analysis of 
the volume data is based on an iterative two-step algo- 
rithm: firstly, the approximate positions of the parti- 
cle centroids are computed through the local maxima of 
cross-correlation with an inscribed sphere with inradius 
Ri — a/(2\/6). Secondly, the orientation and exact po- 
sition of each particle are determined using a steepest 
ascend gradient search. Detected particles are removed 
from the volume and the process is repeated with ran- 
domly chosen starting positions, until no more matches 
are possible. This approach succeeds in detecting more 
than 99.8% of all tetrahedra (final reconstruction illus- 
trated in Fig. Ic). More details can be found in the sup- 
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plementary materials. The further analysis is restricted 
to an inner cylindrical region with diameter 11a and 
height 11a, containing 5056-6051 particles (depending on 
(j)). The vertical variation of the local packing fraction in 
this core region is smaller than 0.002. 

Structural correlations. - A structural measure used 
both in jamming and in liquid theory is the pair correla- 
tion function g{r), which counts the number of neighbors 
at a given distance r, normalized by the binning volume 
and the average density. Fig. 3a shows g{r) for a range 
of (/)■ The first peak is located at r = 0.44a (dashed line) 
and grows monotonely towards denser packings, indicat- 
ing an increase of face-to-face contacts. 

The peak location is close but not identical to the min- 
imal possible distance of 2Ri = rmin ~ 0.408a. Unlike 
the g{r) of sphere packings, the left shoulder of the first 
peak has intrinsically a finite slope. This feature is a con- 
sequence of the lower probability of perfect face-to-face 
alignment compared to slightly shifted or tilted face-to- 
face contacts. One consequence of this specific shape of 
the peak is that compressed packings of soft tetrahedra 
will not show the same square root scaling of contact 
number with volume fraction as observed in sphere pack- 
ings [2. 

g{r) decays quickly after the first peak, showing no 
long-range translational ordering. This is in contrast to 
the g{r) obtained from Monte Carlo (MC) simulations of 
tetrahedral?! [T3] which revealed long-range structure at 
densities as low as = 0.46. The difference can be ex- 
plained by the preparation: the MC simulations prepare 
packings without gravity and friction by slow compres- 
sion of a liquid phase. In contrast our frictional samples 
start to jam due to gravity at (/) « 0.47, after which geo- 
metric frustration and pressure in the pile restrict further 
alignment. 

The orientational order can also be quantified with an 
angular correlation function F(r) (depicted in Fig. 3b I, 
which measures the relative orientation of face normals 
[T5] . F{r) reaches 180° for a perfect face-to-face config- 
uration and 150° for a random configuration. While the 
translational order in g{r) decays quickly, the behavior 
of F{r) suggests long-range orientational order, partic- 
ularly for the densest packings. This result is again in 
disaccord to the findings in [15] but agrees with the nu- 
merical results in [B]. 

Contact numbers. - In order to study the origin of 
mechanical stability, we determine the average contact 
number Z of the tetrahedra. As we can not distinguish 
real contacts from close neighbors based on local force 
measurements |21j . we determine Z from a model based 
fitting method introduced by Aste et al. [55] for sphere 
packings. The supplementary material provides more de- 
tails on this approach which does account for the finite 
resolution of tomography and detection. 

The resulting contact numbers at different densities 
and for the different preparation protocols are shown 
in Fig. [4j The packings prepared by TAP show Z in- 
creasing monotonely with (p. The contact numbers of 
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FIG. 3. (a) Pair correlation function g(r) for different packing 
fractions, (b) Angular face correlation F(r) (Color online) 



the two densest packings, RAMPr=o.5 Etnd VIB, remain 
constant compared to the densest TAP sample. How- 
ever, the packings prepared by RAMP(r = 10) deviate 
distinctly: For comparable 0, contact numbers are ap- 
prox. 0.5 lower, confirming the influence of the prepara- 
tion protocol. Also the initial loose packings show varia- 
tions in Z of the order 0.5. 
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FIG. 4. Average contact number Z as a function of packing 
fraction (p for different preparation protocols. (Color online) 

Constraints. - Contrary to spheres, Z is not suf- 
ficient to evaluate the distance to isostaticity, because 
the number of mechanical constraints fixed by a contact 
depends on the specific contact geometry. Tetrahedra 
have 4 different types of contacts: face-to-face (F2F) con- 
tacts, which are mechanically equivalent to 3 individual 
point contacts, edge-to-face (E2F) contacts (equivalent 
to two point contacts) and the vertex-to-face (V2F) and 
edge-to-edge (E2E) contacts which correspond to a sin- 
gle point contact. For considering the mechanical sta- 
bility, we have to assign one normal force component to 
each point contact. Additionally, friction introduces two 
tangential force components shared by all the in-plane 
point contacts. As these force components are shared be- 
tween two tetrahedra, we obtain the constraint numbers 
Cf2F = 2.5, Ce2F = 2, and Cv2F = Ce2E = 1-5. Mul- 
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tiplying these numbers with the according type-specific 
contact numbers gives the number of constraint per par- 
ticle C. For an isostatic packing C needs to be equal to 
the degrees of freedom each particle possesses, in the case 
of tetrahedra this number is 6 ESI ■ 
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FIG. 5. Top: Number of different contacts per particle as 
a function of packing fraction. The number of F2F and E2F 
contacts increases, while the point contacts (V2F and E2E) 
reach saturation values. Bottom: Variation of the number of 
constraint per particle with packing fraction. (Color online) 

We determine the contact type using the adjacent face- 
to- face-angle aF2F of each pair of contacting tetrahedra 
We e.g. classify a contact as F2F \i apF is in the range 
176.9° to 180°. The lower threshold of this range is based 
on an experimental calibration: We prepare a sample 



with 80 single face-to-face contacts and measure the cu- 
mulative distribution of app in this particular sample. 
This distribution is well described by an error function 
with a point of inflexion at m = 178.2° and a standard 
deviation a ~ 1.3°, suggesting an optimal threshold at 
m — cr. Please see the supplementary material for full 
data and methods. 

Figure [5] shows that the number of F2F and E2F con- 
tacts increases monotonely with 0. At the same time 
the point-like contacts E2E and V2F exhibit a maximum 
around </> = 0.52 or saturate after an initial increase. The 
number of constraints per tetrahedra is between 11 and 
15 and therefore much higher than the isostatic limit of 
6. While this result does not contradict the isostatic 
packings found in simulations of frictionless tetrahedra 
[51 [131 m], there is a clear disagreement with the ex- 
periments using frictional dice reported in [15]. There, 
however, the authors use the contact-specific constraint 
numbers pertinent to frictionless particles. If the con- 
tacts are treated as frictional, the result will also be a 
hyperstatic packing. See also the discussion in [TS]. 

Conclusion. - We have shown that frictional tetrahe- 
dra packings exhibit a hyperstatic constraint number C . 
While the geometry-independent contact number Z sat- 
urates at (f> K, 0.59, the number of constraints C grows 
further, attributable to point contacts being transformed 
into face-to-face and edge-to-face contacts. We find the 
values of C and Z to depend both on packing fraction 
and on preparation history. Our experiments involve 
non-spherical particles and consider friction explicitly, 
thus providing experimental support for extending the 
jamming paradigm. Further work will address differ- 
ent preparation methods to increase the packing fraction, 
and investigate local structures, as F{r) of dense pack- 
ings indicates structures beyond nearest neighbors. 
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